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Abstract: We present a new model for simulating exclusive final states in minimum-bias 
collisions between hadrons. 

In a series of papers we have developed a Monte Carlo model based on Mueller's dipole 
picture of BFKL-evolution, supplemented with non-leading corrections, which has shown to 
be very successful in describing inclusive and semi-inclusive observables in hadron collisions. 
In this paper we present a further extension of this model to also describe exclusive final 
states. This is a highly non-trivial extension, and we have encountered many details that 
influence the description, and for which no guidance from perturbative QCD could be found. 
Hence we have had to make many choices based on semi-classical and phenomenological 
arguments. 

The end result is a new event generator called DIPSY which can be used to simulate complete 
minimum-bias non-diffractive hadronic collision events. Although the description of data 
from the Tevatron and LHC is not quite as good as for PYTHIA8, the most advanced of 
the general purpose event generator programs for these processes, our results are clearly 
competitive, and can be expected to improve with careful tuning. In addition, as our model 
is very different from conventional multiple scattering scenaria, the DIPSY program can be 
used to gain deeper insight in the soft and semi-hard processes involved both in hadronic 
and heavy ion collisions. 
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1. Introduction 

In higfi energy collisions the high density of gluons at small x imply that typical events 
contain several hard partonic subcollisions, and that saturation effects become important. 
A dynamical model based on BFKL evolution [1,2] for small x and including saturation 
effects, has been presented in a series of papers [3-7]. It is based on Mueller's dipole cas- 
cade model [8-10] , which is a formulation of the leading logarithmic (LL) BFKL evolution 
approximation in impact-parameter space. In our model we also include the most essential 
beyond LL corrections to BFKL. In a MC simulation these effects are resummed to all 
orders (see the discussion in ref. [11]). The transverse coordinate space is particularly suit- 
able for treating unitarity constraints, multiple collisions, and saturation effects. Mueller's 
model includes multiple subcollisions in the frame used for the analysis, with similarities 
to saturation effects in the Color Glass Condensate [12-16] or the BK equation [17, 18]. 
In our cascade model we also include saturation effects corresponding to pomeron loops 
within the evolution, in a way which is almost independent of the Lorentz frame used. Also 
effects of confinement are included and important for the result. 

In the dipole picture the colour charge in a gluon is screened by a corresponding 
anticharge in a neighbouring gluon. Gluon radiation implies that one dipole is split in 
two dipoles, and the result of the cascade is a chain of colour connected gluons. When 
two cascades collide, the can interact via gluon exchange. This implies an exchange of 
colour, producing two colour chains forming a BFKL ladder between projectile and target. 
Multiple interaction and saturation represent multiple pomeron exchange and pomeron 
loops. Note here that what is multiple interaction in one Lorentz frame is related to 
saturation within the cascade in another frame. 

The BFKL formalism is directly applicable to inclusive observables, and the cited 
articles include results for total, (quasi)elastic, and diffractive cross sections in pp scattering 
and in DIS, in very good agreement with experimental data. The numerical results are 
obtained with a Monte Carlo (MC) event generator called DIPSY^. 

In the present paper we want to generalize the model to describe exclusive final states. 
Our aim is here not to generate the most precise predictions, and compete with traditional 
generators like PYTHIA8 and HERWIG++. Instead we want to test our understanding of 
small X evolution and saturation effects. In the traditional programs the evolution to small 
X is determined by structure functions, for which the behaviour at a cutoff scale Qq is 
tuned to experimental data, and saturation effects are tuned by the value of a cutoff p±o 

^Acronym for Dipoles in Impact-Parameter Space and rapidity (Y) 
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Figure 1: A chain of initial state gluon emissions. We denote the real emitted ghions in a gluon 
ladder q±i, and the virtual links k±i. 



for hard subcoUisions, determined individually for every energy. In our model these effects 
are instead generated dynamically, based on a perturbative BFKL pomeron and screening 
effects from a dipole picture, including a model for pomeron loops. In addition, in our 
approach we can also study nontrivial effects of fluctuations and correlations, which are 
neglected in most analyses. In the present paper we study pp collisions, but the model 
can also be applied to interactions with nuclei, where it gives a unique possibility to study 
finite size effects and effects of fluctuations. We want to return to these problems in future 
publications. 

In order to generate exclusive states it is necessary to take into account colour coherence 
and angular ordering as well as soft radiation. The latter includes contributions from the 
z = 1 singularity in the gluon splitting function. These effects are taken into account in 
the CCFM formalism [19,20], which also reproduces the BFKL result for the inclusive 
cross section. In ref. [21], describing the Linked Dipole Chain model, it was demonstrated 
that the inclusive cross section is fully determined by a subset of the gluons in the CCFM 
approach, denoted "A;_|_-changing" gluons. As illustrated in fig. |l], we denote the real emitted 
gluons in a gluon ladder q±i, and the virtual links fc^j. Momentum conservation here implies 
that k^j_i = k^j + q_|_j. In a A;_L-changing emission the difference is large between two 
adjacent virtual links, k±i_i and k±i. This means that q±i ~ max(A;_|_j_i, Softer 
emissions with smaller q±i have k±i_i ~ k±i, and are "/c_|_-conserving" . In CCFM gluon 
emission is associated with "non-Sudakov" form factors, and in ref. [21] it was demonstrated 
that these gluons can be summed up in such a way, that the total cross section is fully 
specified by the A;j_-changing emissions with no associated non-Sudakov form factors. The 
remaining A;j_-changing gluons (/^™ are called "primary gluons" in ref. [21] and "backbone 
gluons" in ref. [22]. They are ordered in both lightcone variables g+ and g_, and therefore 
also in rapidity or angle, in accordance with QCD coherence. (In the LDC model [21] 
emissions with k±i ~ on the boundary between A;_|_-conserving and A;_|_-changing 

emissions, were treated slightly different from the BFKL prescription, which also implied a 
slightly different value for the exponent A in the power- like increase oc In the present 

model the LL BFKL result is exactly reproduced.) 
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The chains of primary gluons also determine the structure of the final state, but as 
discussed above, to get the full exclusive states softer emissions must be added as final 
state radiation. In summary the generation of exclusive final states contains the following 
steps: 

1. The generation of two dipole cascades from the projectile and the target, in accor- 
dance with BFKL dynamics and saturation. 

2. Calculating which pairs of one parton in the projectile and one parton in the target 
do interact. In BFKL the emission of gluons is a Poissonian type process, and the 
interaction probability is calculated in the eikonal approximation. 

3. Extracting the primary, A;_|_-changing, gluons, and checking that they are given the 
correct weight. These gluons form colour connected chains between the projectile and 
the target, including branchings and loops. This step includes the removal of branches 
in the cascade which do not interact, and restoring energy-momentum conservation. 

4. Adding final state radiation in relevant parts of phase space. The result consists of 
chains of colour connected gluons. 

5. Hadronization. Here we use Lund string hadronization, where colour strings are 
stretched between the colour-connected gluons. 

The outline of the paper is as follows: In section Q we describe our dipole cascade model 
for inclusive observables. The problems encountered when generating exclusive final states 
are discussed in section ^. In section Q we describe the key procedures used in the MC 
for solving the problems presented in section |3|, while some of the more technical details 
are left for the appendices. Some results and predictions are given in section ^, together 
with a discussion of tunable parameters in the program. We conclude with a summary and 
outlook in section ^. 



2. The Lund dipole cascade model for inclusive cross sections 
2.1 Mueller's cascade model and the eikonal formalism 

Mueller's dipole cascade model [8-10] is a formulation of LL BFKL evolution in transverse 
coordinate space. Gluon radiation from the colour charge in a parent quark or gluon is 
screened by the accompanying anticharge in the colour dipole. This suppresses emissions at 
large transverse separation, which corresponds to the suppression of small kj_ in BFKL. For 
a dipole {x,y) the probability per unit rapidity (Y) for emission of a gluon at transverse 
position z is given by 

dV a 2 {x- yf . , , _ 3a^ 

—— = — d z- TTTT 777, With a = . (2.1 

This emission implies that the dipole is split into two dipoles, which (in the large Nc limit) 
emit new gluons independently, as illustrated in fig. |2[ The result reproduces the BFKL 
evolution, with the number of dipoles growing exponentially with Y . 
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The evolution of the dipole cascade in transverse coordinate space. In each step, a 



Figure 2: 

dipole can split into two new dipoles with decay probability given by eq. (2.1) 




proj. 




targ. 




Figure 3: A dipole-dipole interaction implies exchange of colour and reconnection of the dipole 
chains in the colliding cascades. The arrows indicate the direction of the dipole going from colour 
charge to anticharge. 



In a high energy collision, the dipole cascades in the projectile and the target are 
evolved from their rest frames to the rapidities they will have in the specific Lorentz frame 
chosen for the analysis. Two colliding dipoles interact via gluon exchange, which implies 
a colour connection between the projectile and target remnants, as indicated in fig. ^. In 
the Born approximation the interaction probability between one dipole with coordinates 
{xi^Uj) in the projectile, and one with coordinates (xj^yj) in the target, is given by (the 
factor 2 in the definition of / is a convention) 



^fij = ^f{xi,yi\xj,yA = ■ ' 



log 



{xi-Xj)'^{yi-yj) 



(2.2) 



At high energies the strong increase in the number of dipoles gives a large probability 
for multiple dipole-dipole subcollisions, and the transverse coordinate space is particularly 
suitable for the treatment of multiple interactions and unitarity constraints. Assuming 
that the subcollisions are uncorrelated, multiple collisions are taken into account in the 
eikonal approximation, where the probability for an inelastic interaction is given by 

Int. prob. = 1 - e"^^, with F = Y^ fij. (2.3) 

The multiple interactions produce loops of dipole chains, corresponding to the pomeron 
loops in the reggeon formalism. 

Assuming also that the elastic scattering amplitude, T, is driven by absorption into 
inelastic states, we find via the optical theorem 

T = 1 - e"^, (2.4) 
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and thus the probabihty for an elastic interaction given by T^. (We use a definition such 
that T in this case is purely real.) 

To get the final result for the proton-proton cross sections we have to take averages 
over the projectile and target cascades, and integrate over impact parameter h. For the 
total non-diffractive cross section we get 

^Tinei = / d^h{l - e-2^(^)) = j <fh{l - (1 - T{h)f). (2.5) 

When the projectile has an internal structure, which can be excited, the purely elastic cross 
section is obtained by taking the average of the amplitude T, before taking the square: 

aei = j (fh{T{b))\ (2.6) 

Taking first the square gives the total diffractive scattering (see ref. [23]): 

<Tdiff = j cPh{T{hf). (2.7) 

The cross section for diffractive excitation is consequently given by the difference 

CTdiff ex = <Tdiff - (Tel = j Cpb {{T{bf) - {T{b)f) . (2.8) 

Thus diffractive excitation is determined by the fluctuations in the scattering amplitude. 
It is also possible to calculate e.g. the cross section for single diffractive excitation of the 
(right-moving) projectile by taking the average over the target cascade before, but over the 
projectile cascade after, squaring the amplitude. We also here have to subtract the elastic 
scattering, and thus get 

asD = I d'b{{{T{b))l)R-{T{b)f), (2.9) 

where {■)l(r) refers to the average over the left (right) cascade only. 
2.2 The Lund dipole cascade model 

The Lund dipole cascade introduces a number of corrections to Mueller's original formu- 
lation, which are all beyond the leading logarithmic approximation. The corrections are 
described in greater detail in previous articles [3-5], and a short summary is presented 
here. 

2.2.1 Beyond LL BFKL evolution 

The NLL corrections to BFKL evolution have three major sources [24]: 
Non-singular terms in the splitting function: 

These terms suppress large z-values in the individual parton branchings, and corre- 
spond approximately to a veto for Inl/z < 11/12. It prevents the child from being faster 
than its recoiling parent. Most of this effect is taken care of by including energy-momentum 
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conservation in the evolution, and ordering in the hghtcone momentum p+. To first approx- 
imation this corresponds to a veto for z > 0.5 [21]. This cut is effectively taken into account 
by associating a dipole with transverse size r with a transverse momentum k± = 1/r, and 
demanding conservation of p+ in every step in the evolution. In addition this gives an 
effective cutoff for small dipoles, which eliminates the numerical problems encountered in 
the MC implementation by Mueller and Salam [25]. 

Projectile-target symmetry: 

A parton chain should look the same if generated from the target end as from the 
projectile end. The corresponding corrections are also called energy scale terms, and is 
essentially equivalent to the so called consistency constraint [26]. This effect is taken into 
account by conservation of both positive and negative lightcone momentum components, 
p+ and p- . 

The running coupling: 

This is relatively easily included in a MC simulation process. The scale in the running 
coupling is chosen as the largest transverse momentum in the vertex, in accordance with 
the results in ref. [27]. 

It is well known that the NLL corrections are not sufficient to give a realistic result for 
small X evolution. The characteristic function x{l) has spurious singularities at 7 = and 
1, the point 7 = 1/2 does not correspond to a saddle point, and the power Aefr, determining 
the growth for small x, is negative for = 0.2. As discussed in refs. [11,28,29], these 
problems can be cured by a resummation of subleading logs to all orders, which is executed 
automatically in our MC simulation. 

As shown in ref. [21] the projectile-target symmetry and ordering in p_ shifts the 
spurious singularities in xil) to —uj/2 and 1 -|- a;/2 (when the symmetric scale kik2 is used 
in the Mellin transform). Forshaw et al. have demonstrated that if this shift is taken into 
account, the sensitivity to the exact size of the veto in z is very small [30]. Salam has made 
a detailed study of different resummation schemes, and finds that the result is very stable, 
estimating the uncertainty to about 15% [11,24]. 

2.2.2 Non-linear effects and saturation 

As mentioned above, multiple interactions produce loops of dipole chains corresponding to 
pomeron loops. Mueller's model includes all loops cut in the particular Lorentz frame used 
for the analysis, but not loops contained within the evolution of the individual cascades. 
As for dipole scattering the probability for such loops is given by Og, and therefore for- 
mally colour suppressed compared to dipole splitting, which is proportional to a = NcOs/tt. 
These loops are therefore related to the probability that two dipoles have the same colour. 
Two dipoles with the same colour form a quadrupole field. Such a field may be better 
approximated by two dipoles formed by the closest colour-anticolour charges. This corre- 
sponds to a recoupling of the colour dipole chains. We call this process a dipole "swing" . 

Thus double dipole scattering in one Lorentz frame corresponds to a swing in another 
frame. This is illustrated in fig. ^, which shows the collision of two dipole cascades in 
two different Lorentz frames. The horizontal axis represents rapidity, and the vertical axis 
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rapidity 



Figure 4: Creation of the same loop diagram in different interaction frames, shown as the vertical 
dashed lines. The horizontal axis represents rapidity, and the vertical axis the two transverse 
coordinates. The left side shows a central Lorentz frame where multiple interactions (dipolc (23) 
interacting with (46) and dipole (02) interacting with (57)) forms the loop. On the right side, the 
same loop is formed in an asymmetric frame with a single interaction, but with a swing between 
(34) and (25). 

shows S3mibolically the two transverse coordinates. The left column illustrates the collision 
in a central Lorentz frame, where the rest frame is given by the dashed vertical line, while 
the right column shows the same process in an asymmetric frame. Row (a) shows the 
starting configuration with two dipoles. In the left column row (b) shows the evolved 
dipole cascades just before the collision, row (c) the situation when dipoles (23) and (46) 
have interacted, and row (d) after a simultanous interaction between (02) and (57). In 
the right column row (6) shows how the left cascade has a longer development, while in 
this example there is no emission in the right dipole. It is now assumed that dipoles (34) 
and (25) have the same colour, and row (c) shows the cascade after a swing between these 
dipoles. Finally in row (d) dipoles (01) and (35) have interacted. We see that the final 
result is the same. A double interaction in one frame is equivalent to a swing and a single 
interaction in another frame. We note that in this formulation the number of dipoles is 
not reduced. The saturation effect is obtained because the smaller dipoles have smaller 
cross sections. However, when as in this case, the dipoles in the loop (247) in the right 
figure on row (c) do not interact, the process can also be interpreted as a 2 — 1 process, 
which corresponds to the emission of gluon 2 from the right cascade as seen in the central 
Lorentz frame. 

In the MC, dipoles with the same colour arc allowed to swing with a probability 
proportional to {rfr2)/{r'^r1), where ri and r2 are the sizes of the original dipoles and 
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and r4 are the sizes of the recoupled dipoles. Dipoles with the same colour are allowed 
to swing back and forth, which results in an equilibrium, where the smaller dipoles have 
a larger weight. The given weights favour the formation of smaller dipoles, and although 
this scheme does not give an exactly frame independent result, the MC simulations show 
that it is a fairly good approximation. 

2.2.3 Confinement effects 

Confinement effects are included via an effective gluon mass, which gives an exponential 
suppression for very large dipoles [5]. This prevents the proton from growing too fast in 
transverse size, and is also essential to satisfy Froisart's bound at high energies [31]. 

2.3 Initial dipole configurations 

In DIS an initial photon is split into a qq pair, and for larger the wavefunction for 
a virtual photon can be determined perturbatively. The internal structure of the proton 
is, however, governed by soft QCD, and is not possible to calculate perturbatively. In 
our model it is represented by an equilateral triangle formed by three dipoles, and with a 
radius of 3 GeV~^ ~ 0.6 fm. The model should be used at low x, and when the system is 
evolved over a large rapidity range the observable results depend only weakly on the exact 
configuration of the initial dipoles, or whether the charges are treated as (anti)quarks or 
gluons. 



3. From inclusive to exclusive observables 



In this section we discuss the necessary steps when going from inclusive observables to 
exclusive final states. The technical details of the MC implementation are left to sec. Q 
and the appendices. 



3.1 The chain of fc^-changing gluons 

In the introduction we discussed briefly how the structure of an exclusive final state is 
determined by a backbone of gluons. The first step is therefore to extract this backbone 
chain of fe^-changing or primary gluons, with their proper weights. As discussed above, 
the different subcollisions are assumed to be uncorrelated, and the interaction probability 



for two dipoles, i and j, is given by 1 — exp(— 2/jj), with fij given by eq. (2^). When 
the interacting dipoles are determined, the backbone chains can be found by tracing their 
parents and previous ancestors backwards in rapidity. A resulting chain is shown in fig. ^a, 
and we here use the notation qi for the real emitted gluons, and ki for the virtual links in 
the evolution. The same chain is presented in fig. ||b in a {y,lnq'j_) diagram, with y being 
the rapidity. In this plot lng-|- and lng_ (= y it Infe^) increase towards the upper right 
and upper left corners. Due to the "consistency constraint" [26], which is part of the NLL 
corrections to BFKL, the primary gluons are ordered not only in q^, but also in g_. This 
also implies an ordering in rapidity or angle, in accordance with QCD coherence. The real 
emissions satisfy the relation qi+qi- = qfj_, and are represented by points in fig. ||b. The 
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space-like momenta for the virtual links fcj are not constrained by such a relation, and are 
represented by horizontal lines indicating their transverse momentum. 

The weight for a particular chain with (real) gluons qi is given by [21,22] (assuming 
k±o = k±n = 0) 

2— d2/i()(i;q^j), with a = . (3.1) 

i ^ 111 ^ 

We note in particular that this expression is symmetric under exchange of projectile and 
target. The result can also be expressed in terms of the virtual links ki. We then have 

n n—1 

n (fquK^c^u) = n d^^^i^ (3-2) 

i=l i=l 

and transverse momentum conservation implies that we have two different cases: 

1. Step up: ku ~ 9±i > k^-i, with weight d?qu/(l\i ^ 

2. Step down: k±i <C q_Li ~ with weight (Pq±i/q'j_^ ^ d?k^i/k'\_-_^ = d?k±i/k'\_,- x 

Thus steps down in k± are suppressed by a factor k'\_^/k\^_^. This implies that links with 
a local maximum /c_l, as k^ in fig. ^, is given the weight cPk^^i/k^^^ and can be interpreted 
as a hard scattering between /c2 and A;4, producing the two high-p_L gluons (73 and 54. For a 
link with k± lower than the adjacent real gluons, like kj^ in fig. ||, the transverse momentum 
of the real gluons will be determined by the neighbouring links ^3 and k^. Thus in such 
cases the associated weight will be cPk±i, which is non-singular for small k±. Note that no 
particular rest frame is specified for the collision between the two cascades. The result is 
the same if the "hard subcollision" is part of the projectile or target evolution, and thus 
independent of the Lorentz frame used for the analysis. 

3.2 Reabsorption of virtual emissions 

A dipole cascade generates a chain of dipoles linked together by gluons. A dipole-dipole 
interaction via gluon exchange implies exchange of colour, and the chains are recoupled as 
shown in figs. ^ and § Multiple interactions produce dipole loops, as seen in fig. ^ In 
fig. P we show a more complicated example, in which two cascades are evolved in rapidity 
up to the dashed line in the center. Here three dipole pairs interact, forming two dipole 
loops. In a projectile cascade with a large lightcone momentum p+, and initially small 
P-, the partons need a contribution of p- momentum from the target, in order to come 
on shell. Therefore the branches denoted B and C in fig. |6| have to be treated as virtual, 
and to be reabsorbed. The remaining gluons are coming on shell as real gluons, and are 
colour-connected along the solid lines in fig. |^. A colour loop formed by a swing, as the 
one denoted A, can come on shell also if it does not contain any interacting dipole. Such 
a process produces a closed colour loop not traversing the collision line at y = 0. A more 
detailed example of this effect can be found in appendix |C.2| . The gluons in the virtual 
branches reabsorbed, the momenta for the real backbone gluons must be recalculated. This 
process is discussed in more detail in sec. 
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(a) (b) 

-< Ins *■ 

Figure 5: (a) A parton chain stretched between projectile and target, (b) A backbone of 
/c^-changing gluons in a (y,lnq^) plane. The transverse momentum of the virtual links ki are 
represented by horizontal lines. 

3.3 Giving proper weights to the emissions 

In the cascade the gluons at the ends of a dipole are given opposite transverse momenta 
k± = 1/r (see further sec. 13.4] ). In eq. ( |2.1D an emitted dipole of length r is given a weight 
containing the factor cPr/r'^. However, if this dipole emits further dipoles, the new weight 
is proportional to r^. Thus the associated weight is just (Pr for dipoles which have split 
and been replaced by new dipoles. In the cascades shown in figs. |^ and |9| these dipoles 
are marked by dashed lines. In the following they will be referred to as "inner dipoles". 
The remaining dipoles, formed by colour-connected gluons, are marked by heavy lines, and 
they all get a weight proportional to 1 /r^ . They will be referred to as "outer dipoles" . 

Let us study the cascade shown in fig. |^, obtained after the absorption of the virtual 
gluons as described in the previous subsection. The cascade starts from the dipole (01), 
followed by emission of gluons 2, 3, etc. Here the dipoles are first smaller and smaller, 
a ^ b ^ c ^ d. The corresponding /cj^-values therefore become larger and larger in each 
step. After the minimum dipole, with size d, the subsequent emissions, 5, and 6, give again 
larger dipoles with correspondingly lower k_i_ values. The weight containing factors 
for all "remaining dipoles" is proportional to 

d^r2 (Pr3 d^ cPtq J_ 

In this event gluons 3 and 4 recoil against each other with large transverse momenta 
^max = l/d. As all factors of d have canceled in eq. (^^), this gives the proper weight 
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rapidity 



Figure 6: Collision of two dipole cascades in r-rapidity space. The dashed vertical line symbolizes 
the Lorentz frame in which the collision is viewed. The dipole splitting vertex can result in the 
formation of different dipole branches, and loops are formed due to multiple sub-collisions. The 
loop denoted by A is an effect of saturation within the cascade evolution, which can be formed via 
a dipole swing. Branches which do not interact, like those denoted B and C are to be treated as 
virtual, and to be absorbed. 



1 6 




Figure 7: A dipole cascade in r^-space, in which a chain of smaller and smaller dipoles is followed 
by a set of dipoles with increasing sizes. The initial dipole between points and 1 is marked by 
long dashes. Those dipoles which have split into two new dipoles, and thus disappeared from the 
chain, are marked by short dashes. The shortest dipole (34) corresponds to the maximum k±, and 
represents a hard sub-collision. 

cPr^ cx (i^kmax/^max- This reproduces exactly the proper weight for a hard gluon-gluon 
scattering. 

The same cascade is shown in fig. |8| in the {y, ln(/c^)) diagram of fig. |5|b. In this figure 
new emissions are first denoted by open circles, while filled circles mark their position after 
recoil from a later emission. Thus gluon 1 starts at the lower right corner. After the 
emission of gluon 2, it is recoiling to the position marked 1'. Similarly gluon 2 is shifted 
to position 2', when it emits gluon 3. When gluon 4 emits the softer gluon 5, its recoil is 
small and it keeps its position in the diagram. 

Figure ^ shows instead a chain with increasing dipole sizes up to a maximum value, 
''max) which thus corresponds to a minimum transverse momentum, femin- Here we also 
get the correct weight d^rmax/^'max dP^min, which implies that there is no singularity for 



small /cj^-values. Note that all gluons connected to the long dipoles are also connected to 
shorter dipoles, which determine most of their transverse momentum. Thus no gluon has 
q± represented by A;min- 




Figure 9: A cascade where the dipole sizes increase to a maximum, and then decrease. Here the 
size of the largest dipole, denoted 6, corresponds to the minimum k±^ in the chain. 

3.4 Going from transverse coordinate 
space to momentum space 

Cascade evolution 

The dipole picture is formulated in trans- 
verse coordinate space, but experimental data 
are given in momentum space. We must 
therefore make a translation of our results 
from transverse coordinate to transverse mo- 
mentum space. Note that most of the inclu- 
sive observables studied previously have been 
dominated by large dipoles, while the final- 
state observables are typically dominated by 
the large g^, originating from the small dipoles. Thus effects such as proton size and con- 
finement will be seen to affect observables very little, while new phenomena important for 
small dipoles will come into play. The one important exception to this is the deep inelastic 
cross section for large Q^, which gauges small dipoles with an inclusive observable. 

The weight d^r/r^ for gluon emission in the dipole picture corresponds to d'^k_\_/k\ 
in BFKL [1,2] or DGLAP [32-34] evolution. We also found that a hard scattering, which 
in momentum space is proportional to <f'k^/k'^, in the dipole picture corresponds to d'^r. 
These relations are consistent with the way we associate a dipole r with a transverse 
momentum k = r/r^ (which also implies that r = k/Zc^). Thus, although the fixed relation 
between r and k± is inconsistent with Heisenberg's uncertainty relation, it still reproduces 
the correct evolution also in momentum space. This result would also be obtained by 
the relation k± = c/r, for any constant c. (Any result for a cross section would get the 
same number of extra factors of c in the numerator as in the denominator.) We make 
the choice of c = 1 because it gives the correct (A:^) for a Gaussian distribution in r: 
A two-dimensional distribution tx exp(— r^/i?^) has the average (r^) = i?2 The Fourier 
transform of the amplitude oc exp(— r^/2i?^) is oc exp(— k^i?^/2), and thus the density is 
tx exp(-k2i?2) and (k^) = l/R^. 




Figure 8: The same cascade as in fig. 
The X direction shows rapidity, and y direction 
shows In(fc^). The diagonal line is the constant 
of the incoming particle. 
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Dipole-dipole interaction 

It is, however, not possible to make the same identification in the dipole-dipole inter- 
action cross section from gluon exchange. The logarithmic expression in eq. (p. 2D originates 



from a two-dimensional Fourier transform of the propagator l/k'j_, including interference 
from interactions between the equal and unequal charges in the two dipoles. If we replace 
r by l/k± in eq. (2.2), we would obtain a cross section which grows like {dkl/kl)-ln'^k± 



for large k±. Thus, taking the Fourier transform of the scattering amplitude, and then go 
back to momentum space via the relation k± = 1/r does not give back the original cross 
section in momentum space. It is then more correct to use the original scattering cross 
section directly in momentum space. 

The amplitude for scattering of an elementary charge against a target charge at position 
Ttargct corresponds to the (two-dimensional) Fourier transform of the logarithmic Coulomb 
potential in two dimensions, denoted ^(r — rtarget): 

Aik) = ^Jd^r e^'^"-F(r - rtarget) « • (3.4) 

The phase is here inessential for the cross section, which is determined by The phase 
is, however, important when the target is a dipole. For scattering against a colour charge 
at r3 and an anticharge at r4, we get 

A{k) oc ^{e'''^^ - e*'^'-*) = ^2sm{k^-^^)e'''^^-^+^^^/^ . (3.5) 

2 

For a dipole with charge and anticharge at respectively ri and r2, scattering against a 
dipole with coordinates and r4, we get 

A (X -j^4:sm{pi ■ k)sin{p2 • k), (3.6) 

with pi = (ri — r2)/2 and p2 = (rs — r4)/2. Averaging and summing over colours this 
gives the scattering cross section 

^ = 8^ sm^ip.k) sin^p^k). (3.7) 

To generate k-values according to this distribution we simply select random uniformly 
distributed impact-parameters rjnt and make the identification k = ''"int/^fnt- With d^nnt = 
d'^k/k^ this corresponds to the interaction probability 

fij = 8a; sm^( 2 ) sm ). (3.8) 

Note that for small rjnt, fij goes as a constant times the rapidly oscillating sine func- 
tions, with average (1/2)^. This corresponds to the expected d'^k±/k^ for large k±. In the 
limit where the interaction distance is large compared to the dipole sizes, fij falls off as 
rf^^, which gives the infrared stable result oc d'^k± for small k±. 
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3.5 Final state radiation and hadronization 

As discussed above, the backbone chain of /c_L-changing gkions determines the inclusive 
inelastic cross section. To get the exclusive final states the soft /c^-conserving emissions 
have to be added as final state radiation [19,21,22]. As the A;_L-conserving emissions have 
Q±i < k±i ~ k±i-i, they fill the area below the horizontal lines in fig. |5|. Final state 
radiation is also emitted in the jets produced by the backbone gluons, here represented by 
the folds sticking out of the {y,lnq'j_) plane. The separation of e.g. the colour charge in 
gluon g3 and the corresponding anticharge in gluon q2 forms a colour dipole, and gives a 
gluon cascade similar to the cascade in an e^e~-annihilation event. The only difference is 
that emissions with q± > k±2 are not allowed. 

The final-state radiation is thus added in much the same way as in the so-called LDC 
Monte Carlo [35,36]. In the simulations we use the p± ordered dipole cascade model [37,38] 
as implemented in the ARIADNE program [39]^. In this formalism a gluon emission splits 
a dipole into two, both of which may continue radiating independently. Thus, also after 
final state radiation, the parton state consists of dipole chains formed by colour connected 
gluons. 

In the final step the partons hadronize into final state particles. This is treated by the 
Lund string fragmentation model [40,41] as implemented in the PYTHIA8 program [42,43]. 
Note that the parameters in ARIADNE have been thoroughly tuned together with the 
string fragmentation parameters to agree with LEP data, and we do not change any of 
these parameters here. 



4. Generating the exclusive final states 

The previous section outlined some of the conceptual problems encountered when describ- 
ing exclusive final states in the dipole model, and how they can be solved. This chapter 
will go more into detail on how the Monte Carlo handles different situations, and how the 
pieces of the model are fitted together. 

The first step is to select the interacting dipoles from the virtual cascade, which makes 
it possible to identify the on-shell and virtual gluons. The on-shell gluons will form the 
backbone gluons in the previous section, while the virtual gluons are removed. It will turn 
out that some of the backbone gluons have to be reweighted and some of the q± maxima 



will be removed, as will be described in sec. [4.3| . Then a further check of ordering will 
be made to match phase space with final state radiation, and the colour fiow between the 
remaining backbone gluons will be set based on the virtual cascade. These gluons will then 
undergo final state radiation, and hadronize to produce the final state. 

Some of the more involved details, such as how the allowed phase space is chosen 
for emissions and interactions, and how saturation effects in the cascade and interaction 
complicate the procedure, will be discussed further in the appendices. 



■^We will here use a preliminary reimplementation in THEPEG of the old ARIADNE Fortran code. 
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4.1 Selecting the interactions 

As in our original Monte Carlo [3], the two incoming particles are first identified with a 
dipole state at their initial rapidity. These starting partons are called the valence partons, 
and are then evolved through rapidity to the interaction frame. Here, a relative transverse 
distance between the two incoming particles, the impact-parameter, is selected, and all 
the dipole-dipole scattering probabilities fij can be calculated. Since the individual dipole 
interactions are assumed to be independent, eq. ( ^ ) can be used to find the probability 
for non-diffractive interaction between dipole i and j: 

l_e-2/». =l-exp(-16a2sin2 ( Bl^pBi) sin' f^^)) • (4.1) 

The interaction distance Tint is set to the distance between two randomly selected partons, 
one from each dipole. 

Just like in the cascade, the kinematics of the interaction is calculated. In this case, 
there will be two partons coming in from the left bringing with them but with a deficit 
in p_, and correspondingly two partons from the right bringing but with a deficit in 
p^. An interaction where the partons do not bring enough momentum to set all partons 
on shell is vetoed, that is fij is set to 0. The allowed phase space will be described more 
in detail in appendix p. 

It should be noted here that although the collision of two dipoles formally is a 4-to-4 
reaction, with two incoming partons from each side, one of the partons in each dipole will 
have emitted the other parton. Thus a dipole interaction only connects the end of one 
backbone gluon chain, but the amplitude depends on the last two partons in the dipole 
picture. 

4.2 Identifying the backbone gluons 

Once the interaction dipoles are selected, one can trace the parents of the interacting 
partons back, and identify the backbone gluons, as illustrated in figure |. The rest of the 
partons will not get the necessary p_ from the colliding particle to come on shell, and will 
be reabsorbed as virtual emissions. 

In a dipole picture, each emission is a coherent sum of emissions from the two par- 
tons at the end of the emitting dipole, and a unique parent cannot be determined. In the 
conventional parton cascade the radiation pattern is separated in two components repre- 
senting independent emissions from the two different charges. The coherence effect is then 
approximated by angular ordering, and the recoil is taken by a single parent. For time-like 
cascades as in e~''e~ -annihilation or final state radiation this is not a dramatic effect. For 
the space- like cascades discussed here it is, however more important if both parents can get 
a recoil and be put on shell, or if the recoil effect for one of them may be so weak that this 
gluon has to be regarded as virtual and be reabsorbed. Naturally either alternative can 
only be approximate, and the question which of them gives the best description cannot be 
answered by perturbative QCD. We have therefore implemented both schemes in the event 
generator, and when comparing the results with experimental data we find best agreement 
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when the parents share the recoil. For the results presented below, we have assumed that 
the relative shares are proportional to where ri and r2 are the distances to the two 

parents. 

4.3 Reweighting outer maxima 




Figure 10: A cascade with a maximum in q±^ where the small dipole corresponding to the large 
q±, does not split. The outer q± maximum is shown in impact-parameter space (left) and y-q± 
space (right). 



Recalling sec. 3.3, two examples showed that giving the outer dipoles a weight of 
cPr/r'^, and the inner dipoles a flat weight d?r reproduced the correct weights for maxima 
and minima in q±. In these examples the minimum or maximum dipole were both inner 
dipoles, that is dipoles that emitted further dipoles. If an outer dipole corresponds to 
a local minimum as in fig. 10, it is given a weight proportional to d?r/r'^. As the two 
connected gluons both are given large transverse momenta q± = 1/r, we get a distribution 
proportional to d?q±/q'j_, giving a too strong tail out to large q'x-values. 

To get the correct weights, these small outer dipoles have to be suppressed by an extra 
factor r^/r?., where r< is the size of the small dipole giving the q± maximum, and r> is the 
corresponding maximum in dipole size later in rapidity. We note here in particular, that 
this extra reweighting was not needed in the calculation of inclusive observables presented 
in earlier articles. These observables depend only on interacting dipoles, which are never 
reweighted or absorbed in this way. 

In the DIPSY event generator the reweighting is implemented by finding and reab- 
sorbing some of the outer maxima, so that the correct weight is restored. The inner and 
outer dipoles are only known after the backbone gluons are identified, and therefore the 
reweighting can be performed only after selecting the interactions and the identification of 
the backbone gluons. This procedure is described in more detail in appendix 



4.4 FSR matching and ordering 

The kinematics for the reweighted backbone gluons can be significantly diff'erent from the 
kinematics in the virtual cascade. To exactly match the phase space covered by ARIADNE, 
the q+ and q- ordering is checked again for the backbone gluons before being passed on 
from DIPSY. 

In Ariadne it is then checked that each gluon in the final state cascade is un-ordered 
w.r.t. the backbone gluons, i.e., that the positive (and negative) light-cone momentum of 
an emitted gluon is less than that of (one of) the backbone gluons from which it is emitted. 
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As the backbone gluons are required to be ordered, all of phase space is covered exactly 
once by the backbone gluons with FSR. 

4.5 Colour flow 

The backbone gluons are chosen and corrected looking only at the parent structure, that is, 
each gluon only remembers which two partons emitted it, independently of any subsequent 
dipole swings, as these are not changing the momenta of the backbone gluons. For final 
state radiation and hadronization it is, however, important to keep track of the colour flow 
between the backbone gluons, when the virtual emissions are reabsorbed. 

This is done by going back to the colour flow of the virtual cascade, and absorb one 
virtual gluon at a time, until only the on-shell gluons are left. When a virtual gluon is 
removed, the two neighboring dipoles are combined to a single dipole. For events with no 
swing, this will always return the colour flow as if the virtual emissions never happened. 
For events with dipole swings the colour flow is more complicated, and will be discussed 
in appendix |^. 

The proton remnant: The starting configuration of a proton is in our model represented 
by a triangle of three dipoles, roughly representing the positions of three valence quarks. 
However, three connected dipoles corresponds to the colour flow of three gluons rather than 
three quarks, and the extra charge may be regarded as representing the contribution from 
the gluons present in the proton wavefunction already at low virtuality. 

This extra charge was introduced for inclusive observables, but for exclusive final states 
it has other effects that have to be handled. The three valence partons will in general con- 
tinue down the beam-pipe, and if all three dipoles have interacted, then each of the valence 
partons will have two colour connections to the colliding projectile. This is an overesti- 
mate of the colour flow to the proton remnant, as the extra colour charges representing 
the gluonic component in general do not carry a large fraction of the proton energy. This 
is corrected for in DIPSY by reconnecting on average half of the colour flow to the valence 
partons by a dipole swing, which will move the gluonic content of the proton one "rung" 
down the gluon ladder. 

4.6 Higher order corrections 

The above outlines the main points of how the backbone gluons are generated in DIPSY. 
There are, however, further details which, despite their non-leading nature, have to be 
accounted for. Below follows short summaries of these corrections, while the detailed 
description of the algorithms implemented in DIPSY are left for the appendices. 

4.6.1 Coherence as relaxed ordering 

Previously in this section it was described how the outer maxima were reweighted, 
and absorbed with a certain probability. A q± maximum will through g_ ordering veto 
emissions at small q± close in rapidity, but if the maximum is reabsorbed, the ordered 
phase space will be significantly larger. This is solved by overestimating the ordered phase 
space in the initial generation, making it possible to emit gluons with lower q± in case 
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the high-gx gluon is absorbed. This extra phase space can be regarded as an effect of 
coherence, in the sense that several close-by partons can coherently emit a gluon with a 
longer wavelength, and the gluon can use the combined energy in all the close-by partons. 
This process is described in detail in appendix |B|, where it is also shown how the ordering 
and energy conservation in the interaction is implemented, based on how it is done in the 
cascade. 

4.6.2 Saturation effects 

Previous considerations were all based on a single chain of dipoles, but at a 7 TeV collision 
there are on average more than three subcollisions, making saturation effects very important 
for the final state. In appendix |^ it is shown how the procedures described here has been 
developed to work also in events with both multiple interaction and swings. Multiple 
interactions give rise to splitting chains of backbone gluons while the swing can cause both 
splitting and merging chains, and together they are able to build any diagram of gluon 
chains. 

In this section it is demonstrated how the reweighting of outer q± maxima give the 
correct weight (fq±a'^{qi_)/qj_ in every situation, with the only exception being a special 
configuration in two merging chains, where one of the running couplings can get an incorrect 
scale. This is estimated to be a negligible effect. 

The colour flow in a saturated cascade is more complicated as the dipole swing will 
reconnect the colour flow, simulating soft gluon exchanges between different parts of the 
gluon chains. 

5. Self-consistency and tuning 

There are several details in DIPSY that are not fixed by perturbative calculations, but 
nonetheless have a significant effect on the results. These details will have to be decided by 
other means, or tuned to experimental data. Many of them cannot be directly related to 
a tunable parameter, but are rather choices between different approximations or models. 
Other effects were tested, only to turn out to not give improved results. Thus, rather than 
trying to list which parameters and choices were tuned, we will in this section describe the 
constraints and data, which were used to fix the details of the model. After this we will 
compare the model to the remaining observables. 

One of the most important self-consistency constraint is the frame-independence, that 
is the property that all observables should be the same no matter what collision frame yo 
is used. This is a symmetry that is necessarily present in an all-order calculation, but can 
not be expected to be exactly manifested in a fix-order treatment as the non-leading effects 
enter differently in the interaction and the cascade. Thus, unknown non-leading effect will, 
if possible, be tuned to fit a known all-order property. 

This section will first describe how the frame independence of different inclusive and 
exclusive observables fixes many of the details mention in sec. ^ and then how tuning to 
experimental data fixes the last uncertainties. The section is ended with a comparison of 
DIPSY with experimental data and other event generators for a selection of observables. 
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5.1 Achieving frame independence 

Here the frame dependence of some observables are studied, and many non-leading sub- 
tleties can be fixed by requiring an approximately flat frame dependence. Notice that at 
this stage no comparison is made to experimental data. 




0.1 I 1 ' < 1 1 < 1 1 1 0.2 0.4 0.6 0.8 1 

-8 -6 -4 -2 2 4 6 8 (y„ - Y )/(Y, - Y ) 

(a) ' (6) 

Figure 11: (a) Two colliding virtual cascades in the y-q± representation, generated by the Monte 
Carlo in an off-central frame. The lines shows the colour flow between the partons, the dashed 
vertical line is the collision frame t/o a-nd the circled partons are the starting valence partons. The 
outside triangle show the incoming and p- of the protons, (b) The frame dependence of the 
total pp cross section at y/s = 200 GeV. The dashed curve has full q± ordering in the interaction, 
the full one only requires enough energy to set the interaction on shell. 



(Tpp frame independence: We find that the restrictions in the interaction have to be 
very generous to maintain a constant total pp cross section for i/q close to the valence 
rapidities. A proton evolved over the full range to the other protons rest frame, at the 
Tevatron this is 15 units of rapidity, will have large transverse momenta, while the valence 
partons of the unevolved proton will have very small transverse momenta. Requiring full 
ordering in lightcone momentum in this case will disallow almost all dipole pairs, and the 
interaction probability will be strongly suppressed in off-central frames. This is illustrated 



in fig. 11a with a sample virtual cascade at the Tevatron. To maintain the cross section for 
all frames, a very open ordering in fij is needed, allowing interaction between partons even 
when they are not ordered from the virtual cascades. The minimum amount of ordering 



from appendix B.2 gives a maximum deviation of about 20% at the endpoints as can be 
seen in fig. |Tl|b. Anything stricter results in stronger frame dependence. The required 
ordering is thus set to this minimum 

q+q- > Wq±int (5.1) 

where q+ and q- are the lightcone momenta the incoming partons bring, and q±int is 
the interaction momentum used in fij in eq. (|3.8| ) . This requirement only enforces that the 
incoming particles have sufficient energy to set the interacting partons on shell, without 
any additional ordering. This is significantly worse than previous versions of this model. 
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where deviations between different frames were within a few percent. This is is a sign 
that either the interaction ordering is too strict, blocking the interaction towards the end 
frames, or that the evolution does not grow fast enough, and a cascade evolved over a long 
rapidity range is not as strong as two cascades, each evolved only to the middle. 

The first option is unlikely, as the or- 
dering is already as generous as reasonable 
can be done, and in fact tests with no in- 
teraction ordering at all give no significant 
improvement. Likely, this problem is a re- I 
suit on the ordering in the virtual cascade, 
and how the coherence is used. During 
tuning it has shown that this frame depen- 
dence is quite sensitive to these choices, o 0.2 0.4 o.e o.e 1 
and further tuning could possibly solve - + - 
this problem. For this tune however, the ^.^^^^ p^^^^ dependence of 7*^ at = 
result is the one shown in fig. [l|b. 14 ^nd W = 220, showing three different orderings 

in the virtual cascade. For low tjq the proton is 
evolved a larger part of the rapidity interval, and 
vice versa. 
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ton starts at a transverse momentum of 
about 0.3 GeV, while the two quarks from 

the photon start with a much higher transverse momentum. Here the key observation is 
that it is the (7+ ordering that limits emissions of gluons with larger q±, while the g_ order- 
ing limits the emission of smaller q±. Thus, with the interaction frame close to the photon, 
the cross section will be set by the protons ability to climb up to the higher q±, dominated 
by the ordering, and conversely, with an interaction frame close to the proton, the cross 
section will be determined by the g_ ordering in the photons cascade. 

This provides an opportunity to use 
self consistency to determine which the or- 
dered phase space in the cascade should 
be. In agreement with the arguments in 



appendix [B.l| , it turns out that a fully 
coherent ordering will allow all emissions 
that may be ordered in the backbone 
chain. This is however done by overesti- 
mating emissions of large dipoles in the 
virtual cascade. This is seen in fig. 12 
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where the coherent ordering is tilted with 
a larger cross section when the dipole is 
allowed to evolve. Here the 7* is incom- 
ing from the negative z-axis with a large 
p+, and the proton in incoming form the 
positive z-axis with a large p-. Ignoring 
coherence provides a symmetric frame dependence (although with the same bump as in 



Figure 13: dN/dr] for pp at 900 GeV with de- 
fault DIPSY, and DIPSY with the "single parent" 
setting. The thick line is collided in midframe 
j/o = 0, the thin line in yo — —2.7. 



- 21 - 



pp), but will cut away many emissions that may be ordered in the backbone chain. 

A combination is used, where the q± of the effective parton is used, but the rapidity 
of the single parent. This maintains a symmetric frame dependence while still covering as 
much as possible of the ordered phase space for the backbone gluons. The g_ ordering is 
also very important in deciding the energy dependence for inclusive observables, such as 
the pp total cross section. As will be seen later, it turns out that this combined ordering 
provides an energy dependence that fits very well with experimental data. 

Frame independence of dN/drj: For this observable to be frame independent it is 
needed that any given rapidity produces the same density of charged particles indepen- 
dently of where the interaction frame is. This means that the q± must behave the same 
coming from the cascade or from the interaction, and it should not depend on whether 
the cascade goes on for long, or if it interacts immediately. All this depends heavily on 
the ordered phase space in the virtual cascade and interaction, and in how the backbone 
gluons are handled. In practice, tuning for frame independence of dN/drj is equivalent to 
balancing the q± in the interaction and in the cascade. If the interaction has too much 
transverse activity, then the charge particles will tend to clump up around the interaction 
frame when it is moved. Conversely, if the cascade provides too much q±, then the side 
with the longer cascade will have more charged particles, and the interaction on the other 
side in rapidity will not be enough to balance this. 

With the virtual cascade and fij fixed from the inclusive observables according to 
the above arguments, it is only the handling of the backbone gluons left to tune. The 
problems connected to the reweighting of the outer q± maxima mainly entered through 
the ordered phase space in the virtual cascade. Once the virtual cascade and interactions 
are determined, the ordering of the backbone gluons is fixed by matching to FSR, so few 
choices are left at this point. 

One choice, though, is whether both or just one parent should come on shell, as was 
described in sec. [4.2| . Keeping both parents will give more activity in the cascade as more 
particles will be kept, and more recoils will be done. Thus, the single parent choice will 
have the multiplicity more towards the interaction side, and the double parent choice more 
towards the cascade side, as is seen in fig. 

Neither of them deviate more than 10% at any pseudo-rapidity which is encouraging 
for the model, but neither of the options is much better than the other, so this observable 
cannot fix this choice. It will instead be tuned to the experimental value of this observable, 
which is close to line for double parent. 

5.2 Tuning to experimental data 

Inclusive pp cross sections: In the last published tune of the Monte Carlo [6], there 
were 4 parameters, set to: 

Aqcd = 0.2 GeV iJ^ax = 2.9 GeV-^ Rp = 3.0 GeV-^ Wp = GeV~\ (5.2) 
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being the QCD scale, the confinement range, the proton size, and the fluctuation in the 
proton size. These parameters were tuned to the total and elastic pp cross section as 
function of ^/s, and the elastic differential cross section as function of t. 

In principle the same can be done now, but one of the four parameters now hold a 
more important role than in previous papers. Recall from the discussion in section |3.4| that 
the large dipoles will dominate the inclusive pp cross section, while the small dipoles will 
dominate the exclusive observables. One should note that i?max, Rp and Wp mainly affect 
large dipoles, while Aqcd also affects small dipoles. Thus the inclusive pp cross section 
can be tuned using .Rmax, Rp and Wp without affecting the exclusive observables much. A 



comparison of the default tune with data can be found in fig. 14 




HERAW=220 ■ 
DIPSY Default 



1000 

\'s (GeV) 



10 
(QeV) 



Figure 14: The inclusive data that was tuned to. (Left) The total and elastic cross section as 
function of -^s for pp. Experimental data from [44-52]. (Right) The total 7*p cross section as 
function of for W ~ 220 GeV. Experimental data from [53-55]. 

It should be noted that further fluctuations have been included in the proton wavefunc- 
tion other than the fluctuation in size. First, the partons are allowed to fluctuate in shape, 
making the normally equilateral triangle of valence partons distorted. This has essentially 
the same effect as the fluctuation is size as a flatter triangle will have a smaller cross section. 
Second, the valence partons have a Gaussian smear in rapidity. While they previously were 
all placed at the same rapidity, that of the proton, they are now a bit spread out. This 
has a small effect also for the exclusive observables, as the peak in dN/drj and dpj_/dr] at 
the valence rapidity becomes more smeared out. As in previous publications, tuning shows 
that all these fluctuations have to be small to maintain the ratio between the elastic and 
total cross section. Although the fluctuations are not needed for any observable, the fluc- 
tuations in size, angle and rapidity are set to a Gaussian distribution of widths 0.1 GeV, 
0.1 and 0.1 respectively. 

Inclusive 7*p total cross sections: As mentioned several times above, this is a key 
observable for tuning. Changes that affects the 7*p cross section will also affect exclusive 
observables such as dN/dr], and changes to the virtual cascade that affects the dN/drj 
will also affect the 'y*p cross section. However, changes to how the backbone gluons are 
handled after the virtual cascade is made, will affect only the final state, not the inclusive 
observables. a^*p is allowing us to tune the virtual cascade for small dipoles separately, and 
after that is fixed, details in how the backbone gluons are treated can be tuned separately. 
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This observable is affected, apart from the ordering that was fixed above, by Aqcd 



and c from sec. 3.4, that is the parameter that sets q±_ = c/r. With c = 1 as was motivated 



previously, the model can be made to agree with HERA data with a very reasonable^ Aqcd 



of 0.23 GeV, see fig. 14 



dN/drj: Now the virtual cascade is completely fixed, and all that is left is to fix the 
remaining uncertainties in how the backbone gluons are selected and handled. The re- 
maining details essentially decide how many gluons will come on shell, and how many will 
be reabsorbed. 



The most important choice left to be fixed is the problem in sec. whether one or 
both of a partons parents should come on shell. Not keeping both parents produce about 
20% less charged particles than keeping both, and comparing experimental data, this is a 
bit too low. Thus we will in this tune chose to keep both parents in each emission. 

This fixes all the details of the tune in this paper, and are ready to compare the tuned 
Monte Carlo to experimental data. 

5.3 Comparison with experiments 

We will now compare the results from our program with experimental data on exclusive 
final-state observables, both from the LHC and the Tevatron. 

The DIPSY program is written in C++ using THEPEG [56], which is a toolkit and 
framework for implementing event generation models. This framework is also used by 
HERWIG++ [57] and, more importantly here, a new version of ARIADNE [39], a pre-release 
of which we have been using for the results presented in the following. The Lund string 
fragmentation model of PYTHIA 8 has been interfaced to THEPEG and is also used with 
DIPSY. In this way we produce exclusive hadronic final-state, which can be directly com- 
pared to data. 

We will use the RIVET framework [58] (version 1.5.0) for validating event generators, 
which is also available from within THEPEG. We have selected some representative ob- 
servables from underlying event, and minimum bias studies by the ATLAS experiment at 
center of mass energies of 900 GeV and 7 TeV [59,60]. We also show some comparisons 
with results from the CDF experiment at 1.8 TeV [61]. 

It should be noted that these observables are very difficult to describe also for the state- 
of-the-art general purpose event generators such as PYTHIA 8 [42] and HERWIG++ [57], 
which can be plainly seen from a recent review of event generators [62] and on the web site 



http://mcplots.cern.ch/ where up-to-date comparisons between event generators and 
data are presented. 

The most advanced of the general purpose event generators, PYTHIA 8, is the only 
one that gives a fair description of most minimum-bias and underlying-event data, and in 
the following we will compare our results, not only to data but also to a recent tuning of 
PYTHIA 8 (called 4C in [63]). In fact we will show two versions of this tuning, one which 
only includes non-diffr active events, and one which also includes elastic and diffr active 
scattering. Several experimental measurements include a correction to what is referred 



'For comparison, Aqcd ~ 0.22 GeV in the final state shower in ARIADNE, after a tune to LEP data. 
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to non-single-diffractive events, which in principle should be directly comparable to our 
results. But this correction is very model-dependent (see e.g. [62] for a discussion on this) 
and we therefore give an indication of the size of this correction by including both non- 
diffractive and diffractive results from PYTHIA8. Note that our results should be compared 
to the non-diffr active results from PYTHIA8. 

Also for DIPSY we show two sets of curves. As discussed above, a serious constraint 
for our model is that it should be independent of the Lorentz frame in which we perform 
the collisions. To quantify this independence we have performed simulations not only at 
central rapidity, where most measurements are made, but also in an asymmetric frame'^. In 
the central rapidity frame, the transverse momenta of final-state gluons that are relevant 
for the observables are predominantly given by the dipole-dipole interactions, while in the 
asymmetric case the transverse momenta mainly comes from recoils in the evolution. 

In the following we will only show small fraction of the plots we have produced with 



Rivet. Further plots can be inspected on tittp : //home . thep . lu. se/~leif /DIPSY. html 




Figure 15: The charged multiplicity distribution at 0.9 (left) and 7 TeV (right). The data points 
are the ones given in RIVET version 1.5.0 and are taken from [60] and include only tracks with 
p_L > 500 McV in events with more than 6 charged tracks. The full lines are the DIPSY results, 
the dotted lines are from PYTHIA8 with diffractive and non-diffractive events, the dashed lines 
are PYTHIA8 with only non-diffractive events, and the crosses are from a DIPSY simulation in an 
asymmetric frame. 



Minimum-bias observables: In figures 15-18 we show some standard minimum-bias 



observables as measure by ATLAS at 900 GeV and 7 TeV. A general observation is that 
DIPSY has a slightly too weak energy dependence. Looking, e.g.., in figure ^ we see that 
DIPSY tend to overestimate high multiplicities at low energies and underestimate them 
at high energies. Note, however that the cross section for e.g. multiplicities between 40 



*The rapidity of the collision frame is taken to be 2.7 for 900 GeV, 3.0 for 1800 GeV, and 3.5 for 7 TeV. 
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Charged particle px at 900 GeV, track p± > 500 MeV, for N^h > 6 



Charged particle at /TeV, track p± > 500 MeV, for N^h > 6 
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Pythia8(ND) 
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Figure 16: The transverse momentum distribution of charged particles at 0.9 (left) and 7 TeV 
(right). Data points and lines as in fig. [l5| . 
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Figure 17: The pseudo-rapidity distribution of charged particles at 0.9 (left) and 7 TeV (right). 
Data points and lines as in fig. ^ 



and 50 increases by a factor 40, and a large part of the increase is described by DIPSY (a 
factor 25). The effect is more clearly seen in the pseudo-rapidity dependence of the average 



multiplicity in figure 17, where again the energy dependence in DIPSY is a bit too weak. 



For the transverse momentum distribution in figure [l^, the energy dependence is a 
bit better, but DIPSY tend to overshoot both at small and large transverse momenta, 
while undershooting at medium values. For large transverse momenta, we believe that 
our description may improve if we include 2—7-2 matrix element corrections for local q± 
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Charged (pi> vs. N^h at gooGeV, track px > 500 MeV, for N^h > 1 



Charged (pi> vs. Nj, at /TeV, track px > 500 MeV, for Nd, > 1 
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Figure 18: The average transverse momentum as a function of charged particle multiphcity at 
0.9 (left) and 7 TeV (right). Data points and lines as in fig. |l5[ 



maxima, both in the evolution and in the interactions. For the softer part of the spectrum, 
we believe that the discrepancy is related to the point-like valence partons in the proton 
wavefunction. Even when evolved to mid-rapidity, the valence partons will still have a 
large g+ and will allow emissions at large qj_. Even though these emissions will be in mid- 
rapidity, the balancing recoil will go with the valence partons down the beam-pipe, making 
the spectrum in the main ATLAS detectors too soft. 

In the transverse momentum distribution we also see that DIPSY in the asymmetric 
frame gives a much harder spectrum, indicating that the transverse momentum generated 
in the virtual cascade is harder than the what is generated in the dipole-dipole interactions. 



The frame dependence is also clearly seen in the pseudo-rapidity distribution in figure 17 



where we see an increased activity for large rapidities in the asymmetric case, where the 
particles mainly originate from the evolution rather than the interactions. 



In figure 18 we see that DIPSY underestimate the average transverse momentum at 



small and medium multiplicities, which is related to the underestimate of the number of 
particles with medium transverse momentum. 



Underlying-event observables: In figures we show observables related to the 

underlying event, many of which were introduced by Field [64]. The emphasis is to look at 
particle flow and transverse momenta away from jets and/or high transverse momentum 
charged particles to as much as possible isolate effects of the underlying events. 

In figure ^ we show the multiplicity and summed transverse momentum as a function 
of the azimuth angle from a leading charged track with transverse momentum larger than 
2 GeV. We see that all programs have difficulties reproducing the activity close to the 
trigger particle (the towards region) and in the opposite {away) region which should be 
dominated by the recoil from the leading track. In fact, all programs have similar problems. 
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Figure 19: The multiplicity (top) and scalar sum of transverse momenta (bottom) of charge 
particles as a function of azimuth angle w.r.t. a leading charged particle of at least 2 GeV transverse 
momenta at 0.9 TeV (left), and 7 TeV (right). The data points are the ones given in RIVET version 
1.5.0 and are taken from [59]. The generated data as in fig. nR 



especially in the towards region, and while DIPSY in an asymmetric frame is clearly the 
worst, giving too much back-to-back activity, the default DIPSY actually gives a reasonable 
shape. 

It is, however, the transverse region, defined as 60° < A(j) < 120°, which is the most 
sensitive to the underlying event. In figures 2C and 21 we show the number and summed 
transverse momenta of charged particles in this region as a function of the transverse 
momentum of the leading track. We see that the PYTHIA 8 here gives a better description of 
data, while DIPSY as before seem to has a bit too weak energy dependence and consistently 
seem to overestimate the activity at small transverse momentum of the trigger particle. 
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Transverse N^hg density vs. p^^, \/s = 900 GeV Transverse N^^g density vs. = 7 TeV 
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Figure 20: The multiplicity of charge particles in the transverse region, as a function of the 
transverse momentum of the leading charged particle at 0.9 (top left), 1.8 (bottom) and 7 TeV (top 
right). Data points and lines as in fig. ^ We also show data from CDF [61], but note that this 
data is differently normalized. 



6. Conclusions 

We have here presented a completely new model for exclusive final states in non-diffractive 
hadronic collisions, and implemented it in a Monet Carlo program called DIPSY^. 

As mentioned in the introduction our aim is not to produce the most precise model 
for proton collisions, but rather to study the features of small x evolution and effects of 
saturation. The model is based on BFKL evolution, including resummation of sub-leading 

^The DIPSY program is dependent on the Ariadne program for final-state radiation. While this part 
of Ariadne is well tested and working, other parts of this program is still being worked on and Ariadne is 
therefore not officialy released. The official release of DIPSY will be done together with that of Ariadne. 
Meanwhile pre-releases of both DIPSY and Ariadne are available on request from the authors. 
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Figure 21: The scalar sum of transverse momenta of charge particles in the transverse region, as 
a function of the transverse momentum of the leading charged particle at 0.9 (top left), 1.8 (bottom) 
and 7 TeV (top right). Data points and lines as in fig. |20|. 



logs, and saturation effects which includes pomeron loops not only due to multiple interac- 
tions but also within the cascade evolution. In earlier publications we have demonstrated 
that a single perturbative pomeron can describe not only total and and elastic scattering, 
including the scale dependence of the effective power XesiQ'^) [4], but also diffractive ex- 
citation in pp collisions and DIS [5,7]. In this paper we have extended the model to also 
describe exclusive final states in non-diffractive hadronic collisions. 

The model is dramatically different from conventional multiple-interaction scenaria 
implemented in state-of-the-art general purpose event generators such as PYTHIA8 and 
HERWIG++, which are based on DGLAP evolution with structure functions adjusted to 
experimental data. Also on the technical side, the differences are large. We use a forward 
evolution in impact-parameter space and rapidity, while conventional programs use back- 
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ward evolution in momentum space. Our procedure has the advantage that colour screening 
and saturation effects are easily treated. In addition we have a detailed description of fluc- 
tuations and correlations in the proton wave function, which are typically averaged over 
in conventional approaches. The fluctuations are essential in diffr active excitation, and 
correlations are important e.g. in estimates of the effective cross section in double-parton 
scattering [65]. 

On the down-side, although the evolution of inclusive observables is rather well con- 
strained by the resummation of sub-leading logs, the properties of exclusive states suffer 
from a number of non-leading effects, which influence our results. Small dipoles have little 
influence on inclusive cross sections, but a large effect on the final state. We have here had 
to make several choices, which in many cases lack guidance from perturbative QCD, and 
had to be determined by semi-classical and phenomenological arguments. 

Throughout the the development of our model the requirement that the results should 
be independent of the Lorentz frame in which the interactions are performed has been a 
very severe constraint. In addition, our choices have been constrained by the requirement 
that our description of inclusive and semi-inclusive observables, such as the total and elastic 
cross sections remain satisfactory. 

Although a big part of the framework from previous publications was reevaluated to 
describe exclusive observables, the part of the model that describes the inclusive observables 
remains basically unchanged. The reweighting of the final state gluons does not affect the 
interaction probability, and the ordering in the virtual cascade and interaction, which were 
tuned for frame-independence, are very similar to what we have presented in previous 
publications. 

Thus, while the DIPSY program contains a fair amount of parameters and switches, 
maybe even as many as in the multiple-interaction model in PYTHIA8, once the frame 
independence and inclusive cross sections are accounted for, there is little left that will in- 
fluence the exclusive final states produced. As an example, the dependence on the collision 
energy of final-state observables is locked to the energy dependence of the total and elastic 
cross sections, and thus determined by the BFKL evolution and saturation effects. In the 
case of the multiple-interaction model in PYTHIA8, in contrast, the energy dependence is 
a free parameter, although indirectly related to the increase of the total cross section with 
energy. 

6.1 Outlook 

Although the basic cascade evolution appears to be fairly stable, the generation of exclu- 
sive states introduces larger uncertainties, where we have to be guided by semiclassical 
arguments. Also, as pointed out above, the Lorentz-frame independence requirement is 
not quite satisfied. With this in mind, we are very pleased that although not perfect, our 
model gives quite reasonable results for minimum-bias and underlying-event observables 
at the Tevatron and the LHC. This implies that our model can be used to study if tuned 
parameters in conventional event generators, like e.g. the low-p± cutoff for parton subcol- 
lisions, agree with expectations for low- a; evolution and saturation, and if large effects of 
fluctuations are expected to be important. 
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For future improvements we believe the most important are the matching to matrix 
elements for the hard subcollisions, and the inclusion of quarks. The latter is particularly 
important for EW processes, like Z, W, and Higgs production, and for forward production 
where our simple proton model is not expected to be sufficient. 

We also plan to apply our model to collisions involving nuclei, where saturation effects 
are particularly important. We are here able to produce complete partonic final states in 
eA, pA, and AA collision, including finite size effects and non-trivial correlations between 
partons in different nucleons in the nuclei. Our model does not include any hydro-dynamical 
evolution, but in case a quark-gluon plasma is formed, the parton state can be used as initial 
condition for a hydro expansion. In the absence of plasma, e.g. in eA, pA, or collisions 
between light nuclei, we are also able to generate full hadronic states. Also in collisions 
between heavier nuclei it can be interesting to study the collective effects by comparing 
the results with and without hydro. 
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Appendices 

A. max reweighting in DIPSY 



As was described in sec. 4.3, the maxima in in q± corresponding to an "outer" dipole in the 
cascade, that is a dipole that have emitted no on-shell emissions, come in with the wrong 
weight, and need to be reweighted. 

This is in DIPSY obtained by finding and reabsorbing some of the outer maxima to 
restore the correct weight. The inner and outer dipoles are only known after the backbone 
gluons are identified, and thus the reweighting can be performed only after selecting the 
interactions and the identification of the backbone gluons. 

To find the emissions that need to be reweighted, we note that an outer dipole corre- 



sponds to a parton that has not emitted on one side in colour flow, as parton 3 in fig. 10 
The maximum in q± is associated with the the backwards dipole being smaller than the 
forward dipole, that is, dipole c being smaller than dipole d. 

Starting from the interaction frame, the backbone gluons are checked for partons fitting 
these criteria, and the ones found are reweighted by an extra factor 

« = 4^. (A.1) 

where the running couplings have been included for better accuracy. This changes the 
weight of the two dipoles from VFbefore to Waiter, with 



1 



W^before = d^r<d^r> 'Y' , Waiter = d'r<d'r>^\^ — . (A.2) 
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In momentum space, this corresponds to 

as{q>) as{q<:) 



Wbefore = (fq>(fq<- 



■2„ ,2„ o-'i{q>] 



ql 



q\ 



q> 



(A.3) 



Notice that the smah dipole size r< is associated with the large recoil g>, and conversely 
for the large dipole. The reweighting is done by reabsorbing partons of this type with a 



probability 1 — i?, as is illustrated in fig. 22. This reabsorption is done to the outer parent, 
so that their recoil will be canceled out, and the large q± vanish. 





rapidity 



1 



Figure 22: The backbone gluons after the outer q± maximum in fig. plj|has been merged. The left 
figure is in impact-parameter space. 



Some care has to be taken with the definition of 
an outer q± maximum close to the interacting frame, 
where it is not enough to look at the cascade from just 
one side. A dipole interacts with a probability roughly 
proportional to r^, canceling out the f/r^ weight as- 
sociated with its emission, and the colour flow will be 





Figure 23: Partons f to 3 incoming 
from the left side collides with a cas- 
cade from the right, with partons 4 
and 5. Dipole a interacts with dipole 
/, which makes dipole c an outer q± 
maximum. The figure is in impact- 
parameter space. 



recoupled over to the other state. Looking at fig. 23, 
parton 3 will not fit the above definition for outer qj_, 
as it has no emissions on either side. However, since 
the dipole on one side is interacting, and the dipole on 
the other side is not, dipole c will still be a q± maxi- 
mum with incorrect weight. The Monte Carlo uses a 
generalized algorithm that covers also these cases, and 

uses the dipole d in fig. ^ as the larger dipole to reweight and possible reabsorb parton 3. 

Outer q± maxima from the interaction, for example if d would be smallest dipole, does 
not have to be corrected as the interaction amplitude fij already has the correct weight 
for large q±. 

This procedure ensures that the spectrum of large qj_ from the cascade and interaction 
is indeed going as (g>) _ "Yhe proportionality constant is set from the known amplitude 

for a strictly increasing chain in qj_, and the known "penalty" for being unordered. 

This extra reweighting was not needed before, as only inclusive observables were con- 
sidered which are not affected by this procedure. An interacting dipole is always counted 
as an inner dipole in the final state, which means that no interacting dipole will be removed 
in this way, and the inclusive interaction probability for the two cascades is unchanged. 



B. Absorbed partons and ordering 
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2 



rapidity 



In the previous sections, we describe how the backbone gluons are 
controlled for ordering and reweighted, and some partons are absorbed. 
An absorbed parton returns the g+ to its parents, and the transverse 
recoils are undone. This will in general move the parents back in 
rapidity and down in g_L, which opens up a larger ordered phase space 
for the parents as in fig. 24, After parton 2 is emitted, taking part 
of the g+ from the parent and giving a recoil, 1' is still allowed to 
continue emitting. For these emissions to be ordered in and g_, 
they have to be between the dotted lines in fig. 24 . If parton 2 was not 
emitted though, the larger order phase space between the dashed lines 
would have been ordered with respect to 1. 

If the virtual cascade would only allow emissions ordered in (7+ 
and g_, every absorbed gluon would open a part of phase space for 
emissions from the parents that has not been considered in the virtual 
cascade. This is solved by partially relaxing the ordering in the virtual 
cascade, and thus overestimating the ordered phase space to cover the case of reabsorbed 
partons. Care has to be taken though, as the virtual cascade, before any reabsorbtions, 
will be used to calculate inclusive observables. As the backbone gluons are not yet decided 
when the cascade is made, it is from this approach not possible to use the correct ordered 
phase space already for inclusive observables, and at the same time cover all necessary 
phase space for all possible sets of backbone gluons. 

However, by estimating what the backbone gluons may look like already in the virtual 
cascade one can find a sufficiently good middle way: not too strict ordering so that an 
important part of the possible final states get cut away, but not too open overestimate so 
that the inclusive cross sections infiates. 



Figure 24: 

An emission lim- 
its the ordered 
phase space of the 
parents. 



B.l Coherence 

Three observation at this point will all point towards the same solution: 

• Only large q^ emissions, like the one in fig. significantly change the parents 
momenta. Thus it is only absorption of a step up in q^ that risks shadowing emissions 
from the parent. 

• The phase space that get shadowed mainly concerns small emissions. That is, the 
issue is when a large q^ parton is emitted and later absorbed, where the parent has 
not been allowed to emit with small q^. 

• Many of the absorbed backbone gluons are absorbed due to an outer q^ maximum. 
This is exactly the scenario described in the previous two points, that is an absorbed 
step up in q^ that shadows a step down in q^. 

In DIPSY, when a dipole is emitted in the virtual cascade, it is not only emitted by the 
parton at the end of the dipole, but coherently by all partons being close to the end of the 
dipole. This means that recoil and ordering is using an effective parton with the summed 
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momenta of all the partons within range, where the range is set by the distance from the 
parent to the emitted gluon. 



; 

; 

/ V 




rapidity 




rapidity 



l\ 

I 

I 




rapidity 



:i 



(a) 



Figure 25: A parton emission (4), being emitted at different distances. The upper figures are in 
impact-parameter space, and the lower figures are in y-q± representation, where the dotted lines 
show the limits of the q+ and q- ordering. The circles in the upper figures show the coherence 
range of the emission. 



This procedure is shown in fig. 25 , where parton 4 is emitted at different distances from 
its parent. In (a), the emission is a smaller dipole than the previous one, and the coherence 
range is smaller than the distance to the previous parton 2, so the recoil and ordering will 
be done by parton 3 only. If the emission is at a larger distance than the previous emission, 
as in (6), then the emission will be done by an effective parton with the merged momenta 
of parton 2 and 3. Notice that the recoil between 2 and 3 will cancel when their momenta 
are added, and the effective parton will have the momentum of the original parton 2, as 
shown in the lower figure (6). Emissions at larger distance will continue to add parton to 
the effective parton, and can allow emissions that are not ordered with the parent itself, as 
in (c). 

An outer q± maximum is always a small dipole followed by a larger dipole, and 
reweighting will sometimes merge the two close by partons to one. The extra ordered 
phase space opened by this merging is exactly the part of phase space used by the coherent 

and|5|(6). This 



emissions at long distance, which can be seen by comparing figures 10 



algorithm thus ensures that an absorbed outer q± maximum can not open up a part of 
phase space that was not considered in the virtual cascade. By solving this issue, we have 
been forced to introduce coherence, where an emitted gluon cannot resolve partons that 
are closer together than its wavelength. 

While this increased ordered phase space will cover all situations in the final state, 
some care has to be taken. First, some of the small dipoles that emit coherently will not 



-35- 



be absorbed, and in that case the virtual cascade has emitted large dipoles in unordered 
regions. These large dipoles have a significant effect on inclusive observables, and the total 
cross sections may get inflated. 

Second, the large unordered dipoles are emitted with the motivation that the previous 
small dipole may get emitted. However, the situation can be that the small dipole will be 
absorbed because the large dipole was emitted, as the large emission makes the previous 
small dipole a q± maximum. Thus coherence is an estimate of how the final state will look 
which is to some extent self-fulfilling. 

With this in mind, full coherence as described above may need to be scaled back to 



give a more accurate description. This was addressed in sec. p.2| where frame independence 
was used in deciding exactly what phase space should should be considered as ordered in 
the virtual cascade. 

B.2 Ordering in the interaction 

Like with the virtual cascade it is impossible to use an ordering that allows the correct part 
of phase space for both inclusive and exclusive observables, but again some observations 
can be made. 

1. The final result should not depend on whether a recoil comes from the cascade or 
the interaction. Thus the ordering in the interaction should allow the same backbone 
gluons as the ordering in the cascade. 

2. It is not possible to use the exact same ordering as in the virtual cascade, as the 
situation is fundamentally different. At a certain rapidity in the virtual cascade, the 
partons have not yet gotten recoils from emissions at later rapidities. Thus, ordering 
in the virtual cascade is checked when only recoils with partons in one direction have 
been made. The same is not true in the interaction, where all recoils in both cascades 
are already made. 



3. Events where an outer q± maximum next to an interaction, as in fig. 23, is absorbed 
will allow more phase space for the interaction in exactly the same way as in the 
cascade. In the same way, coherence can be used to make sure all final states are still 
being considered. 

4. In the virtual cascade, g+ and q- held two fundamentally different roles where q+ was 
the lightcone momentum the particle came in with, and g_ measured the increasing 
virtuality. In the interaction the situation is completely symmetric in q^ and g_, so 
the orderings must be the same. 

5. In the cascade, there is no absolute limit on how small q- can be allowed at a certain 
rapidity, as no momentum, g+ or g_, will be emitted in the limit of a soft gluon. 
However, energy conservation sets a hard limit on how much g+ can be allowed, 
namely the g+ of the parent. 

Putting these observations together, we see that there is a limit for how much phase space 
can be allowed, as there always have to be enough energy for the interaction recoil to 
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come on shell. Notice that this corresponds to both g+ and q- ordering, as the particles 
from both states must bring sufficient energy. As in the previous section, the possibility of 
merged partons forces coherent ordering when the interaction distance is larger than the 
previous dipoles in the cascade. The allowed phase space is illustrated for two different 



impact-parameters in fig. 26. 




rapidity 



(a) 




rapidity 



Figure 26: Ordering in an interaction for two different impact-parameters, where parton 1-3 are 
belonging to one cascade, and 4 and 5 to the other. The upper figures are in impact-parameter space, 
and the circles shows the coherence range in the interaction for parton 3 and 4. The lower figures 
are in y-q± representation, where vertical dashed line is the interaction frame yo, the horizontal full 
line marks the interaction q±, and the diagonal dotted lines marks the g_ and (?+ parton 3 and 4 
respectively brings to the interaction. 



In (a), parton 3 and 4 are close to each other, corresponding to a large interaction q± 
that will recoil the partons past each other as can be seen in the lower figure. As 4' needs 
more Q- than 3 can supply, and conversely parton 3' needs more than 4 supplies, the 
interaction q± can not be set on shell, and fij will be set to 0. In this case, there is no 
outer g_L maximum in either cascade as the maximum is in the interaction, so no parton 
risk being absorbed and coherence will not include any extra partons. 

If the impact-parameter between the two cascades would be such as in (b) where the 
partons 3 and 4 are farther away from each other, the interaction q± is smaller and the 
coherence range larger. As can be seen in the lower figure, the smaller virtuality and the 
coherent interaction of partons 2 and 3 means that there is enough g+ and g_ to set the 
interaction on shell, and fij is calculated as normal. 

Here it should also be taken into account that by transferring a fraction 1 — z of the 
supplied to the other state at fix q± will increase the virtuality by a factor 1/z, and 
conversely for the colliding parton, which gives the limit for allowed interactions: 



q+q- > IGgiint, 



(B.I) 



where g+ and g_ are the supphed hghtcone momenta from the interacting effective partons, 
and g^int is the transverse recoil from the interaction. 

This is a minimum for what has to be cut away in the ordering, and from frame 
independence it turns out that it can not be stricter (see section 5.1), thus motivating us 
to set the ahowed phase space in fij to eq. ( p.l[ ) . 



C. Saturation effects 



All motivations and examples discussed above have been based on a single backbone gluon 
chain, but in LHC events there will be several chains branching of and merging due to 
multiple interactions and the swing. This complicates the situation for most of the argu- 
ments presented previously in this section. However, the algorithms have been designed 
with saturation effects in mind, to give realistic results in all situations. Here follows some 
of the larger effects from saturation, and how the above algorithms handle them. 



C.l Multiple interaction 




In the interaction, multiple scatterings are 
allowed, which means that there will be sev- 
eral chains cutting the interaction frame. These 
chains will merge at some point when following 
them back, if not before, they will always meet at 
the valence partons. This does not introduce any 
complications for the ordering. Each emission is 
controlled at the time of the emission, and if an- 
other branch has taken energy from the parents, 
then the ordered phase space for emission will be 
smaller, which is expected. 

One has to be a bit careful with the reweight- 
ing however, but the algorithm described in 
sec. 4^ is made to handle both splitting and merging chains. The dipoles created with a 
weight of cPr/r'^ are still the partons that have on-shell children on one side, but not on 
the other, and therefore the definition of an outer q± maximum can remain. 

This is illustrated in fig. where a backbone chain splits up in two due to multiple 
interaction. Dipoles (24) and (34) both get the weight d?r corresponding to d'^qj_/qj_, and 
they will get two factors of as{r) in case of q± maximum. Thus they do not need any 
extra suppression, and they do indeed not fit the definition of outer dipole, as parton 4 
has children on both sides. For all other partons and dipoles, the arguments for the single 
chain can be applied as normal, so the reweighting algorithm works also with multiple 
interactions. 



Figure 27: Two interactions makes a 
chain split in two. Dotted lines show par- 
ent structure, full lines show colour flow. 
The picture is in impact-parameter space. 



C.2 The Swing 
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The dipole swing in the virtual cascade does not 
cause any recoils, but it does affect the parent struc- 
ture of future emissions. In a non-swinged, single 
interaction chain, one of the two parents of each par- 
ton will also be parent to the other parent, making 
for the typical ladder-like structure seen in previous 
figures. However, with swinged emissions this is no 
longer true, and the two parents can have 4 different 
parents. 

This is seen in fig. ^ where a swing has recou- 
ped dipole (45) and (26) into (24) and (56), (56) 
then evolves further and interacts, while (24) does 
not interact. Here the parents of parton 8 and 9 
cover all the partons (notice that the parents of 7 
are parton 5 and 6), so all partons are set on-shell. 
Notice in the figure that even though parton 4 is 




Figure 28: A swing between (45) 
and (26) causing two chains of back- 
bone gluons to merge. Dotted lines 
show parent structure, full lines show 
colour flow. The picture is in impact- 
parameter space. 



colour connected to 2 from the swing, parton 4 was 

not emitted by 2 (it was emitted by and 1); thus no recoil has been made between the 
two partons, and no reweighting has to be made. 

Dipole (14) has the weight from the cascade and is an outer dipole as parton 4 
has no children on that side. Thus (14) will, if smaller than (45), be reweighted by 



(C.l) 



from eq. (AT), and reweighting again produces the correct weights. Correspondingly for 
parton 2, which is in the same situation as 4, but since 2 is a valence parton with no 
parents, there is no reweighting to be done. Dipole (45) is not an outer dipole because 
parton 7 is a child of 5 on the same side as (45), so no reweighting will be done. If (45) is 
a q± maximum, it will be emitted with a weight a(?'45)/r|5, and the swing is made with a 
giving the weight 



probability ''f^l'^ 



''45a(^45) 



a(?-45) 



(C.2) 



Notice that while the power of qj_ is the correct one, a power in the running coupling is 
missing. This power of as will instead come from the smallest dipole in the emission of 
parton 7 from 5 and 6. Thus, a saturated cascade can give rise to q± maxima where one 
of the powers in Og is at the wrong scale. 

This flaw can possibly be mended by introducing a running Os in the swing amplitude, 
but in the current version of DIPSY it is not implemented. As it is a next to leading order 
factor in the weight occurring only in special configurations of a saturated cascade, it is 
estimated to be a small effect. 



In conclusion, the algorithm in section O will still provide a satisfactory treatment of 
reweighting outer q± maxima. 
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C.3 colour flow in saturated cascades 



For a saturated cascade, the colour flow does not necessarily return to the original structure 
when the virtual emissions are removed according to The swings are changing the 
colour structure, and after all the virtual emissions are reabsorbed, the colour flow does 
not necessarily reproduce the structure without the virtual emissions. An example of such 



a virtual cascade and backbone gluons is shown in fig. 29 




Figure 29: A cascade with swings and multiple interactions. The long dashed lines on the left 
side are the three valence dipoles of a proton. Full lines mark colour flow and dotted lines mark 
emissions. There are three interactions, which are the colour lines passing over the interaction 
rapidity marked by a vertical dashed line to the right. The big dots are the backbone gluons, 
selected through tracing each interaction back towards the valence partons, and the rest of the 
partons are virtual. The left figure is the virtual cascade, and four swings A to D are highlighted 
with thick lines. The right figure shows the backbone gluons after the virtual partons have been 
absorbed. 



Swings where one of the swinged dipoles go on to interact will merge two backbone 
chains, as was seen in sec. p.2| , and are very important for a frame-independent description. 
Swings by an outer dipole, such as A, will not influence the choice of backbone gluons, but 
can reconnect the colour flow between the backbone gluons. A colour reconnection between 
backbone gluons can also happen by a swing between virtual dipoles, if the dipoles originate 
from different parts of the backbone chain, as is the case for swing C. If the two swinging 
virtual dipoles originate from the same part on the backbone chain, as swing D, no difference 
will be seen in the colour flow between the backbone gluons. 

In this way, not only exchange of individual soft gluons between the chains are simu- 
lated, but virtual chains of soft and hard gluons can carry colour flow between the chains, 
or between different parts of the same chain. The virtual cascade can thus give a very dy- 
namic and detailed description on how the colour flow changes in saturated environments. 
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which will be increasingly important at high energies, and can be an essential ingredient 
in heavy ion collisions. 
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